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MATHEMATICS.—On the empirical representation of certain pro- 
duction curves.1. C. E. VAN OrsTRAND, Geological Survey. 


INTRODUCTION 


Production data representing the yearly output from individual 
mines, or from groups of mines included within areas of the order of 
magnitude of states or nations, can generally be represented by a 
curve (a, figure 1) which has a zero value at the origin of time—the 
curve rises irregularly and oftentimes quite abruptly to a maximum 
value and then declines rather slowly; presumably to a zero value, 
in an interval of time which may be assumed to be finite, or infinite. 
The problem to be solved consists of two distinct parts—(a) the selec- 
tion of an empirical equation, and (b) the evaluation of the constants 
of the empirical equation selected to represent the observations. It 
is proposed to treat both of these problems with sufficient detail and 
generality to enable the investigator to adopt at the outset certain well 
defined methods of procedure that will enable him to obtain satisfac- 
tory results with a minimum expenditure of time and energy. 


SELECTION OF AN EMPIRICAL EQUATION 
Expressed in terms of mathematics, we may say that the equation, 
y = FG) 


which represents the production (y) in terms of the time (x) has a 
zero value at the origin and at infinity—between these two values of 
the abscissae, the function is finite, positive, and continuous. The 
first derivative vanishes at infinity, but it may, or may not, vanish 
at the origin. 
1 Published with the permission of the Director of the U. S. Geological Survey. 
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Functions satisfying these conditions can be defined in accordance 
with the theorems of the modern theory of functions. The following 
classification is based on the treatise by Forsyth. 
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Equation (8) consists of polynomial and transcendental expressions. 
The functions, f(x) and P(z, ), are polynomials, but g(x) and gi(-) 


may be either polynomial or transcendental. The exponents in f(z) 
are unrestricted—they may be positive or negative, integral or frac- 
tional. The order of contact at the vanishing points is finite or infin- 
ite according as the singularity is accidental or essential. The 
inequalities tabulated in the column “Limiting Conditions” are 
satisfied when the first derivative vanishes at the origin. Equations 
V, VIII, and IX of the following list are very special cases, 
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*A.R. Forsyru. Theory of functions of a complex variable. Cambridge University 
Press, 1900. A. L. Day andC. E. Van Orstranp. The black body and the measurement 
of extreme temperatures. Astrophysical Jour. 19: 24-30. 1904. 
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respectively, of equations 6, 7, and 4 of the preceding list. With the 
exception of equation VII, the remaining equations in the list are 
easily obtained in a similar manner by introducing zeros at the points 
c, and ¢ instead of the points 0 and ©. Both zeros of equation VII 
are at an infinite distance from the origin. 

Apart from certain relations between the constants, equations I to 
VII, inclusive, represent Pearson’s seven types of frequency distribu- 
tion which have been used with great success by himself and others,’ 
chiefly in the field of biometrics. If we properly define the constants 
in IX and put m = n = 2, the equation represents the distribution of 
molecular velocities in a gas. Putting m = 0, n = 2, the equation 
reduces to Pearson’s type VII, which is the well known normal, Lapla- 
cean, or Gaussian distribution of frequencies on which the classical 
theory of probabilities is based. Pearson’s type III and V are like- 
wise special cases of IX. Two other very special cases of IX, 


y = ax” andy = ae” 


are used to represent decline or extinction curves. These relations 


3W. Pauin Etperton. Frequency-curves and correlation. C. and E. Layton, Lon- 
don. 
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are easily represented by straight lines on logarithmic paper.* Equa- 
tion VIII contains as a special case, the equation deduced by Planck 
to represent at a particular temperature the distribution of energy in 
the spectrum of a black body. 

Pearson has developed criteria based on the method of moments by 
means of which one may select from his list of seven functions the 
particular equation of the group which will most accurately represent 
the observed values. Unfortunately, this method of procedure, at 
least in its present state of development, is not easily adapted either 
to the determination of criteria or to the evaluation of the constants in 
the selected equation when the observed frequency curve is incomplete. 
In view of these defects and the fact that theoretical considerations 
give us practically no information as to the nature of the required 
function, it becomes necessary to adopt a different method of pro- 
cedure. Evidently, variation, or degree of fiexibility, or generality 
of functions can be obtained in only two ways—either by increasing 
the number of constants, or by putting the necessary mechanism into 
the function itself. An extreme case of the first is a Fourier series— 
of the second, a function of functions such as the expressions obtained 
for the representation of a frequency distribution by the method of 
semi-invariants. From the preceding brief discussion, it appears that 
equation IX be.t meets the requirement of possessing a very consider- 
able range of generality expressed in terms of a minimum number of 
constants. It has therefore been selected for a first trial. 


COMPARISON OF METHODS OF EVALUATION OF CONSTANTS OF FREQUENCY 
DISTRIBUTIONS 


The precise evaluation of the constants of a frequency distribution 
oftentimes presents more difficulties than does the selection of an 
equation which is assumed to represent the data of observation. Four 
methods of evaluating the constants are now in use: (1) Method of 
least squares; (2) Pearson’s method of moments;* (3) Thiele’s* method 


*Caartes 8. Larxey. Mathematical determination of production decline curves. 
Mining and Metallurgy 4: 341-344. 1923. Harry M. Rozser. Determining the con- 
stants of oil-production decline curves. Trans. Am. Inst. Mining and Metallurgical 
Engineers. (Birmingham meeting, October, 1924.) Wittarp W. Curier, Jr. Es- 
timation of underground oil reserves by oil-well production curves. Bull. 228 U. S. Bureau 
of Mines, VII + 114 pages, 1924. 

5 Arne Fisuer. An elementary treatise on frequency curves. The Macmillan Co. 
1922, ARNE FisHer. The mathematical theory of probabilities. Vol. 1. Second edi- 
tion. 1923. The Macmillan Co. 




















JAN. 19, 1925 VAN ORSTRAND: PRODUCTION CURVES 23 


of semi-invariants; and (4), Krichewsky’s* method of evaluating the 
constants of the first derivative of the required integral equation. 

Columns 3 and 4 of table 1 give Elderton’s’ computed values 
(ye), and residuals (v) for the data tabulated in columns 1 and 2. 
The values were obtained by adjusting Pearson’s type I by the method 
of moments. Columns 6 and 8 give the computed values obtained by 
adjusting equations IX and I, respectively, by the method of least 
squares. The computed values in column 6 are too large over that 
portion of the curve which approaches the age axis as an asymptote, 
but over the remaining portion of the curve, the residuals obtained 
by this method are less than those obtained by the method of moments. 
Using the arbitrary values, c, = 9.1, c. = 65.0, and adjusting the 
remaining constants of equation I by the method of least squares, we 
obtain the numbers tabulated in column 8. Judged by the usual 
criterion that the sum of the squares of the residuals (=v?) must be a 
minimum, the values tabulated in table 1, namely, 298, 550, and 169, 
show that the method of least squares gives a better adjustment than 
does the method of moments; furthermore, equation IX gives a gradua- 
tion only slightly inferior to that of equation I. 

Column 11 gives the values calculated by Fisher* from the 








equation, 
y = 1130 [¢o(x) + 0.0258 goiti(x) + 0.0158 goir(x)] 
wherein 
1 fz + 0.0195\2 
‘eis = 1 ” 3 ( 1.6240 
? 1.624 V2x ° 


and the remaining terms are respectively, the third and fourth deriva- 
tives of ¢,(x). The constants were evaluated by the method of semi- 
invariants from the data of columns 9 and 10. In column 13 are 
tabulated the values obtained from a least square adjustment of equa- 
tion I, using the arbitrary values, c: = c: = 31.0. The values of m 
and m, were assumed in this case to be independent. Comparison of 
the corresponding values of =v, 124 and 639, suggests that it is prob- 
ably impossible to adjust equation I, either by the method of least 


*S. Kricuewsxy. A method of curve fitting. Ministry of Public Works, Egypt. 
Physical Department Paper No. 8. Government Press, Cairo, 1922, E. W. Wootarp. 
On Krichewsky’s method of fitting frequency. curves. Monthly Weather Review 52, no. 2: 
91-94. 1924. 

7W. Pain Exvperton. Loc. cit. pp. 54-58. 

8A. Fisoer. Loc. cit. First reference, pp. 48-52. 
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squares or by the method of moments so as to obtain the remarkable 
degree of precision obtained in this particular instance by the method 
of semi-invariants. 


TABLE 1.—Comparison OF Metruop or LEAST SQUARES WITH METHOD or MOMENTS 
AND METHOD OF SEMI-INVARIANTS 
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The method of least squares frequently fails to give a satisfactory 
adjustment of Pearson’s type V; after passing the mode, the com- 
puted curve generally falls considerably above the observed curve. 
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This defect can be corrected in part by substitution of equation X for 
equation V. 

The preceding examples serve to call attention to some of the difficul- 
ties involved in obtaining the best graduation of a particular set of 
observations; none of the methods is infallible, and ordinarily, the 
best method of meeting the requirements of a particular problem can 
be found only by trial and investigation. 
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Fig. 1. Curves showing the effect of weighting the observation equations. Curvec 
is weighted; curve b is not weighted. 


ADJUSTMENT OF SOME PETROLEUM PRODUCTION DATA 


The annual production of petroleum® for the years 1875-1923, 
inclusive, for the state of West Virginia is shown graphically in fig- 
ure 1. The numerical values, expressed in terms of millions of barrels, 
are tabulated in column 2 of table 2. 

In figure 1 is shown an anomaly of the least square adjustment. 
Curve 6b was obtained by adjusting equation IX, n = 1, on the basis 
of equal weights of the observation equations written in the form 


log.a + log.« m — xb = log.y. 


®G, B. Ricnarpson. Petroleum in 1919-1921. Mineral Resources of the United 
States, 1921. Part II, pp. 253-333, May 26, 1923. r 
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Curve c shows the result of adjusting the same equation after assign- 
ing to each observation equation the theoretically correct weight,'° 
y*. Following are the values of the constants: 


Curveb. Weight1. log = — 1.4008, m = 1.2776,b = — 0.0376. 
Curvec. Weight y*. log. = —5.6961,m = 4.5267,b = + 0.1609. 


The first value of b is negative instead of positive as it should be. 
Attention is called to this peculiar result as it occurs quite frequently 
in attempting to adjust the complete data of production curves and 
may be incorrectly attributed to an error in the computation. The 
defect disappears if the origin is displaced sufficiently to the right as 
the following results, obtained by displacing the origin to the year 
1889, readily show. 

Weight 1. log.a = 0.19062, m = 1.45976, b = 0.10375, =v? = 73.41. 
Weight y*. log. = 0.80217, m = 1.14872, b = 0.08857, =v? = 65.61. 
Comparison of the values of =v* shows that the best results are again 
obtained by properly weighting the observation equations. Equal 
weights may be assigned to observation equations written in the 
form 


y 2 + y log. 2 dm — xy db = dy = 0 A a <P ane ae 3 (1) 


Using the preceding values of the constants obtained by assigning the 
weight, y?, and solving for da, dm, and db on the basis of equality of 
weights, we obtain the following results: 


Weight 1. log.a = 0.76175, m = 1.16496, b = 0.08935, =v? = 65.06. 


The value of 2v? = 65.06 is practically the same as the former value 
65.61; consequently, the two methods lead to identical results provided 
the values of the constants used in the computation of dy = v in the 
differential equation are a sufficiently close approximation to the cor- 
rect values. 

The true value of the exponent, n, in equation IX is not easily 
determined. In order to obtain distinct types of equations which are 
easily capable of interpretation in terms of the theory of generalized 
probability, Pearson adopted the integral values, + 1, — 1, and + 2. 
Fractional values of n must be admitted, however, if one desires to 


100, M. Leuanpv. Practical least squares. First edition. p. 136, McGraw-Hiil 
Book Co, 1921. 
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minimize =v*. Curves corresponding to various values of n are shown 
in figure 2. All negative values, including the value n = — 1 for 
which the curve is shown in the figure, must be excluded for the reason 
that in all of these cases the first derivative vanishes at the origin. 
Comparison of the values of =v for the remaining values of n shows 
that the true value falls between n = 0 2nd n = 0.5, Values of the 
order of magnitude of n = 0.01 leads to indeterminate solutions for 
the remaining constants in the equation. This result is to be expected 
for in these special cases the value of the exponential factor (e~°") 
tends to approach a constant or a straight line. In the following final 
computations, the value n = 1 has been selected as a sufficiently close 
approximation to the true value. 


TABLE 3.—Constants ror TABLE 2 





Ix xX It! 





(2) (3) (4) 


329.3 329.3 
337.5 337.0 336.7 
57.11 56.86 





1.46607 2.59166 
0.81369 — 
0.07470 0.07470 


0.81373 0.76915 
on 15. 48968 
12.0 12.0 

ae 222.0 
Origin 1890 1890 1902 
Mode 11.65 10.89 —0.93 
Mode—year 1901.65 1900.89 1901.07 

















The constants of the equation, 
m bz 
y =are 


are easily determined by writing the observation equations in the form 


log.a + log. xm — xb = logy. 


As previously stated, theoretical considerations show that each obser- 
vation equation must be multiplied by the corresponding value of 
yo. Using the data contained in the first two columns of Table 2, 
we obtain 33 weighted observation equations of the following form: 
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2.4 log.a + 0.000 m — 2.46 = 2.101 
3.8 log.a + 2.634 m — 7.66 = 5.073 

6.4 log.a + 22.378 m — 211.26 = 11.880, 


the solution of which leads to the results summarized in the second 
column of table 3. The first derivative vanishes at the points 


x= andz = *. 


The latter point is designated mode by Pearson, and peak by the petro- 
leum geologist. 

The ordinates in the preceding computation were placed at the ends 
of the year intervals. Precision is generally supposed to be gained by 
placing the ordinates at the mid-points instead of the end points, 
The results of making the computation in this way are summarized 
in table 2, columns 5, 6, and 7; and in table 3, column 3. In this 
particular case, the value of 2v* is reduced from 59.26 to 57.11, but it 
must not be inferred that a reduction of this character will always be 
obtained, for the value of =v? is dependent also upon the point selected 
for origin. 

Pearson places the initial origin at any convenient point, and then 
solves for ¢c:, the point at which the curve touches the z axis. The 
corresponding procedure by the method of least squares consists in 
solving a differential equation similar to (1) for the differential ele- 
ments da, db, dc, and dm. I have not carried out the computations 
for this solution. The results obtained by transferring the origin to 
the point, x = 12.0, and adjusting type III on the basis of unrestricted 
constants are tabulated in columns 9 and 10 of table 2. Comparison 
of the values summarized in columns 3 and 4 of table 3 shows that the 
two methods lead to identical results. The numerical values are 
represented graphically by the curve designated, type IX, in figure 3. 

Two attempts have been made to apply Pearson’s method of 
moments: first, to the incomplete data (y.) tabulated in table 2, 
column 2; and, second, to assumed complete data obtained by using 
the observed values (y,) (column 2, table 2) from 1890 to 1923, inclu- 
sive, and the computed values represented by that portion of the type 
IX curve, figure 3, which extends from 1923 to the end of the curve. 
The values of « in the respective cases are x = -— 0.0123 andx« = 
— 0.2240, thus indicating that type I should be used in preference 
to type III. Both types fail for incomplete data. Evaluation of the 
constants from the so called complete data leads to the following 
results: 
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- 0.4543 2 12. 1008 
Typel. y = 13.36 (1 ¥ zisu) (1 % new) 


Mode-year = 1889.5 + 10.016 = 1889.516 Skewness = 0.838 
Mean-year = 1889.5 + 24.323 = 1913.823 =v = 63.18 
Standard deviation = 17.07 Zy- = 328.1 


) 1.4749 — 0.0022z 


x 


Mode-year = 1889.5 + 13.473 = 1902.973 Skewness = 0.636 
Mean-year = 1889.5 + 24.323 = 1913.823 yw = 79.22 
Standard deviation 17.07 Lye = 323.3 


Type I appears to be preferable to type III. As the labor involved 
in actually deciding the case would be very great, I have evaluated, 
instead, the constants of type I, using the method of least squares. 
Three computations were made, assuming for c, the arbitrary value, 
12.0, and for c, the successive values, 132.0, 222.0, and 312.0. The 
respective values of =v*, wherein v was carried to tenths of a unit, were 
found to be 60.93, 58.41, and 57.97. The curve which represents 
=v? expressed as a function of c, varies so slowly that these values alone 
do not lead to an accurate determination of the desired constant. 
However, the value c. = 222.0 appears to be a close approximation 
to the true value, and as it does not differ greatly from the value 
2 = 216.4 calculated by the method of moments, it has been used in 
extending the computed curve until it reaches the x axis. The exten- 
sion is designated type I in figure 3. Over the remaining portion of 
the curve, the three types of formulae lead to almost identical results, 
as may be readily verified by comparison of columns 6, 9, and 11 of 
table 2. The values of 2v? tabulated in table 3 differ slightly, how- 
ever, and suggest that types III or [X should be used in preference to 
type I, but the agreement in the values is so close, and departures 
from rigorous methods of evaluation of the constants have been so 
marked, that definite conclusions can not be drawn except to say that 
each of the formulas in question gives a very satisfactory representa- 
tion of the data of observation. 


I! 


SOME GENERAL RESULTS AND SUGGESTIONS 


Table 4 summarizes the results of adjusting the petroleum produc- 
tion data given by Richardson for three states and one western 
field. The constants of the equation 
m — bz 


| 
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were evaluated by the method of least squares from weighted observa- 
tion equations of the general form 
ylog.a + ylog.x«m — ry.b = ylog-Yo. 


The resulting values of the constants are tabulated in columns 1, 
2, and 3. Column 4 contains values of r, the probable error of an 
observation of weight unity. The origins were determined by inspec- 
tion. A correct determination of these points would reduce slightly 


TABLE 4.—Summary or ComputTATIONS FOR STATES AND FIELDS 
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Kern River, 
a ae 2.42293/0. 50828/0. 10208/1.34/1900; 4.98/251.5/241.3/0.042) 63.18/312.4/20.2 




















the values of r. The mode-year is obtained by addition of the modeto 
the year selected as origin—thus, the mode-year or peak-year for the 
state of Ohio is 1884 + 12.26 = 1896.26. 

Pearson’s method of moments equalizes the observed and computed 
areas when applied to complete data. The value, ty. = 328.1, 
obtained by the method of moments from the assumed complete data 
of West Virginia is in close agreement with the observed value, Ly 
= 329.3; on the other hand, all of the values resulting from the appli- 
cation of least squares to the particular cases discussed in this paper 
are too large. The magnitude of the discrepancy is shown in column 
9 of table 4. The error, for example, in computing the total produc- 
tion from 1890 to 1923 inclusive, for the state of West Virginia, is 2.3 
per cent of the observed total production (Zy.) for the same time inter- 
val. The corresponding error resulting from an imperfect application 
of the method of moments is only 0.4 per cent. Unfortunately, this 
advantage is offset to a considerable extent by the failure of the method 
to minimize 2v*; the difference between the two values, 63.18 and 57.11, 
is altogether too large. Mathematicians seem not to have settled the 
question as to the precise relation between the two methods of evalua- 
tion of constants. Jackson" has shown that the methods are identical 


11 DuNHAM Jackson. The method of moments. Am. Math. Monthly 30: 307-311. 
1923. 
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for power series, but the relationship existing between the two methods 
when applied to the more complicated functions has not been estab- 
lished. It is obvious from elementary considerations that the con- 
dition of maximum probability and the other condition of equalization 
of observed and computed areas are not necessarily satisfied at the 








PRODUCTION-MILLUIONS OF BARRELS 











Fig. 3. Curves showing final adjustment of types I and IX 


same time. The computations summarized in columns 3 and 8 
of table 1 provide a numerical proof of this proposition. 

The total area of the computed curve is easily obtained from the 
formula 





anve “d Tr (m+ 1) _ m! 
area = Jaz ¢ ct = yet ~ inti 


The results obtained by substitution in this formula are tabulated in 
column 11 of table 4. The values, 565.3, 792.1, 448.0, and 312.7 
were obtained by summation of the computed ordinates (y.). The 
agreement between the two summations is all that could be desired. 
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As the area of the computed curve from the origin to the mode (area 
abca, figure 3) can not be evaluated by integration; it is necessary to 
base the calculation on the computed values of y. The areas tabulated 
in column 10 represent these summations for the integral years and for 
the fractional part of the year interval extending from the end of the 
last integral year to the mode. The ratio of computed peak produc- 
tion to computed total production is tabulated in column 12. 

The preceding results show that the simple equation (2) can be 
used to’represent a very considerable range of production data with a 
degree of accuracy that is quite sufficient for most purposes. The 
constants are easily determined by the method of least squares, and 
the mode is given by the simple ratio, m/b. The computations repre- 
sented by the curves in figure 3 show that equation I also gives a satis- 
factory representation of the data. For other states or areas, some of 
the remaining equations in the list will undoubtedly prove to be 
useful. 

Pearson’s types of frequency distribution promise to be of great 
value in various fields of science. It may be of interest in this con- 
nection to call attention to the possibility that certain sands, for 
example, may be characterized by the mode, mean, and skewness (see 
figure 3) of the curve which represents the distribution of grain diam- 
eters, and that certain properties of the sands, such as permeability, 
may be related to these characteristics. The results recently obtained 
by Wightman, Trivelli, and Sheppard” on the relation of size-fre- 
quency distribution of grains in photographic emulsions to sensito- 
metric characteristics serve to emphasize the importance of this 
suggestion. a 
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PHYSICS.—Note on the theory of optical dispersion. G. Brett, 
Department of Terrestrial Magnetism. (Communicated by L. 
A. BAvER.) 

Soon after the development of the Spdg = nh principle attempts 
have been made by Sommerfeld and Debye! to account for optical 
dispersion by studying the perturbations caused in the electronic 
motion by the incident wave. These attempts, however, were unsuc- 
cessful because the frequency at which the calculation showed the 
reaction to be strongest proved to be the frequency of revolution of 
the electron and not the absorption frequency given by the relation 
AW =h». This difficulty was avoided by Ladenburg and by Laden- 
burg and Reiche?* by disregarding the nature of-the perturbations and 
by postulating directly that an atom which is in a position to absorb 
radiation of frequency »,. contributes to the refractive index of the 
medium in which it is situated a term which would be contributed by 
a linear oscillator of the same frequency ». A modification of this 
theory has been given by Kramers.* 

In these theories the dispersion is discussed only for frequencies 
which do not lie in the absorption band. For this reason the equiva- 
lent classical resonator is not supposed to have any dissipation of 
energy. It is the purpose of the present note to discuss how far it 
is necessary to introduce the dissipation of energy into the theory. 

On the classical theory a frictional term in the equations of motion 
could be interpreted as due to radiation. On the quantum theory 
there is no continuous radiation. Thus a frictional term introduced 
directly into the equations of motion of the classical resonator would 
appear, to a certain extent, illogical, and one is led to look more closely 
into the possibilities of a classical resonator without a frictional term 
We adopt the notation of Van Vleck,‘ and consider an energy level « 
from which it is possible for the atom to go into higher states r by 
absorbing energy and to descend into lower states ¢ by emitting energy 
Taking the orientations of the atoms to be random, the polarization 
for N, electrons in the state s is hypothesized by Kramers to be 


a c Aye A,_s: | 
P.= 32 x* |. Vrs? (Vrs? — v) ze vee (vee — v*) ad 


1 Despre, Sitzungsber. Miinchener Akademie, p. 1, 1915. 

? LADENBURG and Reicue, Naturwissenschaften 27: 584. July 6, 1923. 

*H. A. Kramers, Nature, May 10, 1924, p. 673; August 30, 1924, p. 310. Seevl 
Van Vueck, Phys. Rev. 24: 344. 1924. 
4 J. H. Van Vick, Phys. Rev. 24: 330. 1924. 
5 Van VLECK, loc. cit. Equation 17. 
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where A,—,,d¢ is the a priori probability that in the time dé an atom 
which is in the state r should fall into the state s. Whether this is 
the actual form or only similar to reality the essence of what follows 
remains unaffected. 

Let us consider one particular dispersion term. 


é Ays 


Pu = 35 Yea? (ne? — *) 





N.€ (1) 


This means that under the action of the incident wave having an 
intensity & the N. atoms which are in the state s suffer a total polariza- 
tion P,. owing to their ability to be transferred to the state r. Now 
it is clear that any single atom can not give a purely monochromatic 
expression (1) because it does not remain in the state s indefinitely 
but on account of the A,,, is bound to leaves spontaneously. There- 
fore, even if one should suppose that while an atom is in the state s 
it gives a polarization of the type (1), one must expect an expression 
different from (1) for the observed average polarization. Expression 
(1) is essentially the reaction on the wave due to a classical oscilla- 
tor and for this reason we shall simplify the discussion by considering 
a resonator having the equation of motion® 


+@2=—B8, (2) 
m 


Here E, is the x component of the incident wave, e, m, 2 are respec- 
tively the charge, mass, and resonant frequency of the resonator. 
The incident wave we suppose to be of the type 


E, = E. cos(wot — ) (3) 


Let us suppose that at the time ¢ = ¢,, the atom falls into the state s 
and that at the time ¢ = ?¢, it falls out of it. It is impossible to say 
(with the present state of our knowledge) what the initial conditions 
of the oscillator (2) are at ¢ = ¢,. Lacking better information, how- 
ever, we may take the oscillator to be in a state of rest and in a position 
of equilibrium at the instant of its birth. The solution of (2) with 
these initial conditions is 


6 While writing this note the author was unaware of the derivation given by H. A. 
Lorentz (Proc. Amst. Acad. 8: 591. 1906.) in which formula (10) is derived. Inasmuch, 
however, as the present note considers the question from the point of view of Fourier 
integrals, and shows the manner in which only the part of the polarization synchronous 
with the incident wave survives, it has not been withdrawn from publication. 
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Ex. 
Rot Se tai (425|- cos (¢ — wots) cos 2 (t — t;) — 


= sin (€ — wot) sin 2 (t — t,) + cos (wot — o| 


In the neighbourhood of #) = @ this is very nearly the same as 


e E.. 


Tires os | eos (wot — €-) — cos {2 (t — th) + wots - i] (4’) 


Letting 
te +i ti > T F (4.5) 


and neglecting  — Qin comparison with 2 we have in the neighbour- 
hood of the absorption frequency the following Fourier Integral for x 





ae a ae <_ Whe — oo) Tf _ sin@ — 2)7 


rm (#—ad)J (LL 2@— 0) ~ Fea |e (ot — 
faa) aes] 
= See [ +. 2 sin(wt — 





@®— Wo _ 
e— (w— wo)t,) } dw (5) 
Let us suppose tha‘. tie intervals (t, ¢.) are distributed throughout a 


long interval .7 in such a way that on the average their density is 1. 
The value which concerns us is 


F 
fz (t; th, T) dt, = 


i =0 


ae is to fens - meal Gre —w)T sin — z *" pony a 
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( sin(“ % “«r) 3 sin'(* . “")) tg 7 a és 


@ — Wn w—Q 





Again, to within the same approximation, the wave (2) when acting 
for a time -7, and analyzed into a Fourier Integral is 


= Ex Be (ene = Se enclshes dhl (7) 


Since 7must be thought of as >>7' the band of frequencies repre- 


sin (w -_ wo) TF 


sented by is to be thought of as infinitely more sharp 


@o-— Wo 


sin(w — wo)T 
than that represented by an ano and consequently we may re- 
Wo 


write (6) as 
< 
f xt tT) dt, = 


“n = 0 


satan) Cea] G - BE aR) ota - 9 - 


sin “5 °r) 
@o—- 
BO ella. ) 
- sin 2 (2 sin(wo — 2)T 


w@ — wo | Re 2(wo — Q) ) sim(at ~ 9 + 


. (wo — 2 
sn 2. ") cos(wt — a| dw (8) 


wo — 
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dz = 0 we may disregard the part of (8) containing 





Comparing what is left of (8) with (7) it is clear that 


@o-— Wo 


the atom acts as though during the time 7 it had a polarization 
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Bal pe ater ee 


e {%o —2 ) 
2 sin{ 2 T 


@o- 


) Cos(wol — «) — 


sin (wot — iY (9) 





If now we should suppose that the times 7’ have a distribution function 


. (i) a(F 


of the type e 
(9) viz. 


7) then we obtain an average value of P from 


€ Bao Te (Q — wo) [To cos(wot — €) + (2 — wo) sin(wot — «€)] 


P be m (Q + wo) {1 a T? (wo — 2)?] 





(10) 


This result for the polarization is thus obtained on the hypothesis of 
intermittent actions, the number of intervals of lengths between T 
iY 


-_ Te T 
and 7’ + dT’ beginning per second being e ™d T, The total number 
9 


of intervals started per second is therefore 1 and the average duration 
of an interval is 7). Hence, if we are to compare (10) with the result 
of the classical theory for a single continuously acting resonator having 
a frictional term then we must divide (10) by 7>. To within the 
approximations made (on performing an elementary calculation for the 
classical resonator) the result is seen to be identical on the two theories. 

We have thus shown that the phase change introduced by the clas- 
sical frictional force may be also thought of quantitatively as due to 
the interrupted action of the resonator. If the Correspondence Princi- 
ple may be applied in this case and if it is to be interpreted as assign- 
ing a certain similarity to the actual dispersion and the theoretical 
classical dispersion then we are practically forced to abandon any other 
influence than that of interrupted actions in the higher quantum states. 
If, further, the width of emission and absorption lines is related to the 
dispersion by the laws of the classical theory then we must think of 
this width as being entirely due to the finite time of each independent 
emission. 

Absorptions beginning in the lowest energy level occupy a singular 
position in this theory because atoms may remain in these indefinitely. 
Whether they suffer phase changes in their motion while in the lowest 
state is unknown but for the sake of uniformity such changes might be 
hypothesized. 














JAN. 19, 1925 DYAR AND SHANNON: NEW MOSQUITOES FROM BRAZIL 39 


It is of interest to observe (as will be brought out in more detail 
in another note) that the above theory of interruptions accounts for 
the polarization of resonance radiation in weak magnetic fields in a 
more satisfactory manner than the classical resonator.’ 


ENTOMOLOGY.—New mosquitoes from Brazil (Diptera, Culicidae). 
Harrison G. Dyar and Raymonp C. SHAnnon, U. S. National 
Museum. (Communicated by 8S. A. RoHWER.) 


Dr. J. Bequaert has placed before us for determination a collection 
of mosquitoes recently made in Brazil, chiefly in the State of Ama- 
zonas. Among a number of interesting species collected, the following 
are hitherto undescribed. 


Sabethoides glaucodaemon, new species. 


Of usual size in the genus, largely purplish black; proboscis extending well 
beyond the antennae, slightly longer than the abdomen; palpi small, slightly 
longer than the two basal flagellar joints; eyes contiguous on lower side of 
head for a greater distance than they are above; prothoracic lobes contiguous 
above, their scales overlapping; rather numerous setae on anterior margin; 
mesonotum with setae only above roots of wings and a few on anterior mar- 
gin; spiracular sclerite with three setae; propleura with two setae; sterno- 
pleura without setae; a small but dense tuft of long setae on upper posterior 
corner of mesepimeron; pleurae with dense white scales below; trochanters 
and base of femora yellow, the under sides of the femora white scaled basally 
along the entire length of the posterior pair; abdomen compressed laterally, 
dark scaled on upper half, yellowish white scaled below, the colors divided 
in a straight line, though the white is illy contrasted in certain lights. Wings 
normal, basal cross-vein opposite the anterior cross-vein; roots of halteres 
yellow, stem and knob blackish; mid tarsi white scaled below on last four 
joints except narrowly at base of second. 

Type, female, no. 27744, U. 8. Nat. Mus.; near San Alberto, Rio Branco, 
Amazonas, Brazil, August 28, 1924 (J. Bequaert). 


Nearest related to imperfectus B.-W. & B.,differing chiefly in the slightly 
longer proboscis, and from both this and chloropterus Humb. in the coloration 
of the abdomen. 


Culex (Choeroporpa) bequaerti, new species 


Rather small dark brown species; occiput with erect forked scales, all the 
recumbent scales broad, mostly white in front, a patch of black ones on each 
side of the middle; antennae fairly long, exceeding length of proboscis, which 
is about equal to the length of abdomen; integument of mesonotum very 
dark brown, scutellum somewhat paler; scales narrow, dark brown; dorsal 
setae sparse but well developed; a row of pronotal setae; pleurae yellowish, 


7In such a way we may think of the field caused by virtual oscillators as virtual in 
the sense that it carries energy away only statistically but does not necessarily produce 
damping of the motions induced in the virtual oscillator by the exciting field. 
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prealar setae about seven, a row of fairly strong setae along posterior margin 
of sternopleura, a little weaker below; a single lower mesepimeral seta; legs 
entirely dark, except for paler ventral surfaces of femora; abdomen dark 
brown, with whitish scales on ventral surface at bases of segments; wing 
scales broad subcostally, those on base of fork of second vein narrowly ovate 
to ligulate; halteres pale, the knobs dark. 

Male: palpi exceeding the proboscis by nearly the length of the last two 
joints, the penultimate joint with a small whitish ring at base, otherwise 
dark. Scales of mesonotum deep bronzy brown; three pronotal setae, broken 
in the female type. 

Hypopygium.—Side piece a little longer than hemispherical; inner division 
of lobe strong, running far into the side piece, with an infuscated patch 
basally, columnar, long, exceeding the outer lobe, with two strong long 
hooked and distorted filaments at tip, one inserted basally of the other; outer 
division small, with four stout rod-like filaments on the oblique outer aspect 
and a small rounded leaf basally of them. Clasper slenderly snout-shaped, 
the spine appendiculate. Tenth sternites comb-shaped, with six teeth, en- 
larged at base, with only a rudiment of basal projection; first plate of meso- 
some normal, the articulated plate rather narrow, emarginate on one side; 
second plate curved, tip furecate, the arms short, inner pointed, outer smooth, 
a long strong horn a little beyond the middle of the stem; basal hooks slender, 
strongly recurved, not projecting at base; ninth tergites conically pointed, 
small, setose, connected by a chitinous band. 

Type, male, paratype, female, no. 27745, U. 8. Nat. Mus.; male, Sororoca, 
Rio Branco, Amazonas, Brazil, September 1, 1924; female, Carmo, with other 
data the same (J. Bequaert). 


The preceding description was written before we had an opportunity of 


examining Miss A. M. Evans’ recent paper' describing new Choeroporpa from 
Brazil. Of the species there described C. (C.) thomasi Evans comes nearest 
to the present form. The mesosomal plate in the two is much the same. 
The comb of the tenth sternite of thomasi appears abnormal in Miss Evans’ 
figure; in bequaerti it consists of seven long equal teeth. The inner division 
of the lobe of the side piece has a longer stem in bequaerti than in thomast, 
the distance between the insertion of the two filaments being less than the 
remaining basal part of the stem, whereas in the figure of thomasi the reverse 
is the case. The outer division of the lobe of the side piece is differently 
formed, having no inner limb in bequaerti and the leaf is inserted on the stem 
basally of the other filaments, whereas in thomasi it arises between the limb 
and the outer setal group. 


Culex (Choeroporpa) albinensis Bonne-Wepster & Bonne 


Culex (Choeroporpa) albinensis Boune-Wepster & Bonne, Ins. Ins. Mens. 7: 
173. 1920. 

Culex (Choeroporpa) albinensis Dyar, Ins. Ins. Mens. 8:62. 1920. 

Culex (Choeroporpa) gordoni Evans, Ann. Trop. Med. & Par. 18:369. 1924. 


No difference is apparent between Miss Evans’ detailed figures of the 


1 Ann. Trop. Med. and Par. 18: 363-375. 1924. 
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hypopygium of gordoni and the characters of albinensis in a slide before us. 
The distribution of the species includes Paramaribo, Surinam, and Mandos, 
Brazil. 


Mansonia indubitans, new species 


Basal antennal joints as dark as rest of antennae; antenna somewhat 
shorter than length of proboscis; palpi as long as four basal flagellar joints; 
proboscis on basal two-fifths with pale and dark scales intermixed, a 
rather broad white ring a little beyond middle, beyond blackish, paler at 
apex; mesonotum dark brown, setae normal, sparse small narrow golden 
scales intermixed, with long dark scales on the sides posteriorly; a number of 
pronotal setae; postspiracular setae present; sterno-pleural on posterior mar- 
gin with one long stout seta midway and smaller setae on either side, mes- 
epimeron with three setae near anterior margin; femora and tibiae with dark 
and light scales intermixed; hind tibiae darker; first tarsal joint without basal 
or apical ring, but with scattered white scales on inner surface; second and 
third tarsal joints white basally, also the fourth joint of mid and hind tarsi, 
remainder dark; wings dark scaled with numerous white ones intermixed, 
all broad; abdomen dark, with triangular patch of white scales on first seg- 
ment, apex directed forward, the venter with numerous broad white scales 
intermixed. Knobs of halteres dark brown. 

Type, female, paratypes, three females, no. 27746, U. 8. Nat. Mus.; Belem, 
Para, September 19, 1924; Carmo, Rio Branco, Amazonas, August 31, 1924; 
above Santarem, July 22, 1924, Itacoatiara, Amazon River, September 15, 
1924 (J. Bequaert). 


Similar to étillans Walker, the palpi shorter, and with slight differences in 
coloration as indicated above. 


Anopheles celidopus, new species 


Medium size, grayish in general appearance; occiput with erect truncate 
white scales above, dark brown below, white setae and scales between eyes; 
antenna shorter than palpus, scales only on the basal flagellar joint; palpus 
but little shorter than proboscis, with outstanding dark brown scales, a few 
white ones at apices of second and third joints; prothoracic lobes with tuft of 
scales above; mesonotum with pale curved hairs sparsely distributed, a little 
denser anteriorly and darker on the sides; pleura with two indefinite pdle 
pollinose lines; legs dark with narrow white rings at apices and bases of all 
but the last tarsal joints; abdomen dark, with sparse dark hairs, a few white 
scales on the dorsum of the last segment, and many dark and white scales on 
the venter of this segment; cerci densely scaled, mostly dark; wings with eight 
more or less definite white spots on anterior margin, and numerous other 
small white spots irregularly distributed over the wing; wing scales lanceo- 
late; knobs of halteres dark. 

Type, female, paratypes, four females, no. 27747, U.S. Nat. Mus.; Carmo, 
Rio Branco, Amazonas, Brazil, September 1, 1924 (J. Bequaert). 


This species does not fit well into any of the existing groups of Anopheles. 
It comes nearest to Arribalzagia, but lacks the lateral scale tufts of the 
abdomen. 
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SCIENTIFIC NOTES AND NEWS 


Dr. S. G. Mortey, Associate in Middle American Archaeology, Carnegie 
Institution of Washington, gave public lectures at the Institution on October 
28 and December 3, on Initial excavations at Chichen Iiza, Yucatan. 
Under agreement with the government of Mexico, the Institution has the 
privilege of excavation and archaeological study at this ancient Maya city 
for a period of ten years, beginning January, 1924. 


Dr. A. E,. Dovetass, Director of the Steward Observatory, University of 
Arizona, gave an illustrated lecture at the Carnegie Institution of Washington 
on December 5, on Cycles in development of tree rings as an evidence of 
climatic variation. 


Two lectures on atmospheric dynamics were given at the Carnegie Insti- 
tution on January 6 and 8, by Pror. V. B)erknes, professor of hydrody. 
nami¢s at the Geophysical Institute at Bergen, Norway, and research 
associate of the Institution. The first lecture was on The forces which lift 
airplanes, and the second, on Problems in dynamical meteorology. 


Dr. J. A. ANDERSON, of the Mt. Wilson Observatory of the Carnegie Insti- 
tution, gave an illustrated lecture at the auditorium of the Institution on 
January 5, under the title: An experimental method of studying high 
temperatures. 











